Granular Brownian Motor. 
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An asymmetric object, undergoing dissipative collisions with surrounding particles, acquires a 
nonzero average velocity. The latter is calculated analytically by an expansion of the Boltzmann 
equation and the result is compared with Monte Carlo simulations. 



I. INTRODUCTION 



Brownian motors are spatially asymmetric construc- 
tions that, operating under nonequilibrium conditions, 
can rectify thermal fluctuations. They have been the 
object of intense study over the past 15 years In 
this letter, we make the connection with another active 
field of research, namely granular matter Q. In such 
systems, particles undergo dissipative colllisions and are, 
by construction, in nonequilibrium. We therefore generi- 
cally expect that a Brownian motor will arise if we break 
spatial symmetry in granular matter. To investigate 
this question in more detail, we propose here a minimal 
model, which has the advantage that the resulting sys- 
tematic speed can be calculated exactly. We bypass the 
difficulties associated to the granular gas itself by assum- 
ing that the collisions between the gasparticles are elastic 
and that the gas is extremely diluted (ideal gas limit). 
Both spatial asymmetry and dissipation are introduced 
by considering an asymmetric object which undergoes 
dissipative collisions with the surrounding particles. The 
outcome of our calculations is that the predicted speeds 
are quite large, namely comparable to the thermal speed, 
and should therefore be easily observable in experiment. 



II. THE MODEL 



An asymmetric object, for example a triangle, is free to 
move, however without rotation, along a horizontal axis, 
see Fig. Q] Its motion, with speed denoted by V, is in- 
duced by dissipative collisions with surrounding particles. 
For simplicity, we consider a two-dimensional system and 
assume that these particles are an ideal gas initially at 
equilibrium at temperature T in an infinitely large con- 
tainer. Hence the spatial distribution of the particles is 
uniform outside the object and their velocity distribu- 
tion is Maxwellian. Note that post-collisional particles 
are no longer at equilibrium. However, in the limit of 
infinite dilution, recollisions with the (convex) object do 
not take place and the hypothesis of molecular chaos (no 
correlations between the speed of the object and that 
of the particles prior to collisions) is valid. The follow- 
ing Boltzmann-Master equation is therefore a microscop- 
ically exact starting point, describing the time evolution 




FIG. 1: Sketch of the system under consideration: an asym- 
metric object moving freely along the horizontal axis as a 
result of the inelastic collisions with the surrounding gas par- 
ticles. 



of the probabilitiy density P(V,t): 

/+oo 
du[W{V - u;u)P(V - u,t) 
-oo 

-W(V;u)P(V,t)]. (1) 

W(V; u) is the transition probability per unit time that 
the object (mass M) changes speed from V to V + u, 
due to collisions with the surrounding particles (mass 
to). Its explicit expression has to be derived from the 
laws of dissipative collisions, taking into account the ge- 
ometrical configuration of the object. To remain con- 
sistent with the molecular chaos assumption, we restrict 
ourselves to convex objects. The latter can be fully char- 
acterized by the shape probability function F(0), defined 
so that F(9)d9 is the fraction of the outer surface with 
polar angle between 9 and 9 + d9 (see Fig. [3]). The to- 
tal circumference will be denoted by S. The object is 
constrained to move along a specific direction, say the 
x-axis, and its velocity is then V = Ve x . The veloc- 
ity components of the particle are denoted by v x and v v , 
and primes will be used to denote post-collisional speeds. 
Since there are no external forces in the x-direction, the 
corresponding component of the total momentum is con- 
served under collision: 



MY = mv x + MV 



(2) 



Furthermore, the interaction force is assumed to be or- 
thogonal to the surface at the point of impact, implying 
that the tangential speed component of the impinging 
particle is also conserved: 



t = v-t 



(3) 



Finally, the dissipative collision reduces the relative ve- 
locity orthogonal to the surface by a factor r (0 < r < 1), 
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FIG. 2: Schematic representation of a convex object. Its 
shape is characterized by the shape probability function F(9). 



The tangent and orthogonal unit vectors at the sur- 
face with polar angle 8 are t = (cos 9, sin 9) and h — 
(sin (2, — cos 9) respectively. The above equations deter- 
mine uniquely the post-collisional speeds in terms of the 
pre-collisional ones. In particular, the change of speed u 
of the object reads: 



u = (1 + r) 



m sin 



m sin z 9 + M 



(v x — V — Vy cot 8) (5) 



the so-called normal restitution coefficient: 
(V' -if)-h= -r(V -v)-h. 



The transition rate VF(F; u), which is the probability per 
(4) unit time for such a change of speed, is thus given by: 



/>2tt r+oo />+oo 

W(V; u)= d9SF(9) / dv x 



dvyO 



(V — v) ■ h (V — v) • n 



x p4>{v x ,v y )5 



u - (1 + r) 



to sin 9 
to sin 2 9 + M 



(v x — V — v y cot 9) 



(6) 



where Q [x] denotes the Heaviside function, p is the density of the gas and 4>(v x ,v y ) the velocity distribution of the 
granular gas, taken to be a Maxwellian: 



4>{vx,v y ) — 2^j, ex P y ^—j 

The integrals over the speed of the colliding gas particles can be performed explicitly, so that VT(V; u) reads: 

r 2TT\ 



W(V;u) = Sp 



2nkT 



pit r \ 

iQ[-u] J +u<d[u] J J A9F{9) 



(to si^^ + M) 2 

x 5— exp 

(1 + r) 2 TO 2 sin 9 

I 



to sin 2 8 ( i u(m sin 2 8 + M) 



2kT 



(1 + r)TOsin 



(7) 



(8) 



The remaining integral depends on the shape probability the calculation of the resulting average drift velocity, we 



function of the object. 



III. RESULTS 



The model introduced above can be considered as a 
granular variant of a thermal Brownian motor Q. For 



can proceed along lines similar to those described in [6j, 
yielding an expansion for the stationary average drift ve- 
locity in terms of the small ratio e — \fmjM. The main 
steps of the calculation are presented in the appendix. 
The resulting expressions for (V) and (V 2 ) read: 
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and 



(v 2 ) = 



+ rkT 



1 kT m , 

1 -r) 

16 M M ' 



(sin 



_4(l + r )2_ 



sin 



(10) 



The brackets involving the polar angle are averages 
with respect to the shape probability function F(6). In 
case of elastic collisions, r = 1, the object reaches a state 
of thermal equilibrium with the surrounding gas, that is 
(V) = and (V 2 ) = kT/M. For dissipative collisions, 
the object obviously does not equilibrate. To lowest or- 
der in e, one verifies that the velocity of the object is 
still Maxwellian, but at an (lower) effective temperature 
Toff — (1 + r )T/2- As a result, there is a continuous flow 
of energy (heat) from the gas into the object. This en- 
ergy flow, in conjuction with the spatial asymmetry, is 
the driving mechanims for its systematic motion. The 
direction of the motion is determined by the geometry. 
Turning to a specific example and for later comparison 



with the simulations, we consider an isocleses triangle 
with apex angle 29q. The shape probability function is: 

2(1 + smffo) V 

+ 5[6-(n-e )}) (11) 

and: 

( S ^0) = tl)^±^. (12) 
x ' l + sin0 o 

The expressions for (V) and (V 2 ) now become: 



(V) 



M 4 V M 2 



/ m V 



7r(sin0o — 1) 



M 2 32 

2(5 + r)(sin0 o - 1)(1 - sin0 o + sin 2 O ) + < 



sin 4 6*0 — 1" 
sin On + 1 



and 



{V 2 ) = i±^— + — — — (1 - r ) [Trfsinflo - l) 2 - 4(1 +r)(l - sin0 o + sin 2 O )1 



(13) 



(14) 



We conclude that the triangle has an average negative 
speed: it always moves opposite to the direction in which 
it points (0o € [0j 7F /2])! The average speed goes to zero, 
consistent with the fact that the motion is produced by 
fluctuations, when the mass ratio m/M goes to zero. It 
becomes a maximum in the limit of high dissipation r — > 
and maximum asymmetry |(sin 3 0)/(sin 2 0)| = 1. 

For comparison with these theoretical results, we have 
performed stochastic simulations to generate trajectories 
of the Master equation. The highly efficient algorithm is 
based upon exact acceptance-rejection methods for gen- 
erating the Maxwellian inflow distribution [7|; see the 
appendix C in [6j for a complete description of the al- 
gorithm. Averages were taken over 40 000 realizations, 
with each realization taking 100 000 time steps. The re- 
sults are shown in Fig. [3J Agreement between theory and 
simulations is excellent. Deviations are observed, as ex- 
pected, when the mass ratio m/M becomes smaller, but 
a surprisingly good agreement persists even for M of the 
order of m, as shown in Fig. 2] 



IV. DISCUSSION 

According to our theoretical result, the speed of the 
granular motor is of the order of the thermal speed 
y/kT /M and should therefore be easily observable. 
There are however a number of critical comments to be 
made concerning its experimental observation. First, one 
needs to address the fluctuations in velocity. Considering 
for simplicity the case of extreme dissipation and asym- 
metry, one finds (V 2 )/(V) 2 = 32/tt x M/m + .... We 
conclude that even for a mass ratio m/M = .1, r.m.s. 
speed fluctuations are still 10 times larger than the av- 
erage speed. Only for extremely light motors (m/M = 5 
according to the now less reliable prediction of the per- 
turbative result) will fluctuations be of the same order 
as the average speed. This difficulty can of course be 
circumvented, either by taking enough sample trajecto- 
ries, or by measuring the displacement for a time much 
longer than the correlation time of the velocity fluctua- 
tions. Second, the average speed is the result of a subtle 
unbalance between the collisions on different sides of the 
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the case of three dimensions. 
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FIG. 4: Comparison between theory (full line) and simulation 
(dots) for the average velocity (V) of the triangle as a function 
of its mass M. Parameter values are set to m = 1, T = 1.0, 
do = 0.2tt, p = 0.00222 and r = 0.95. 
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FIG. 3: Comparison between theory (full line) and simulation 
(dots) for the average velocity (V) of the triangle in a) as a 
function of the angle 8o (T — 1.0) and b) as a function of 
the temperature (8q = 0.27r). In both cases we set m = 1, 
M = 100, p = 0.00222 and r = 0.95. 



object. We therefore expect that the resulting systematic 
speed will depend sensitively on the velocity distribution 
that is used to describe the surrounding gas particles. For 
example, in the case of a driven granular gas, the veloc- 
ity distribution can deviate strongly from the Maxwellian 
form and moreover depends significantly on the thermo- 
stat being used (see [1| for a recent review) . Furthermore, 
a temperature anisotropy can appear in a granular gas 
because of the unidirectional energy input [4j . As a con- 
sequence, experimental observed velocities, using granu- 
lar gases with a non-Maxwellian velocity distribution, are 
expected to yield significant deviations from the theoret- 
ical values predicted here. Deviations from a Maxwellian 
distribtuion can however be addressed by repeating the 
above calculations with the appropriate velocity distri- 
bution. Such calculations reveal that the average speed 
is typically of the same order of magnitude, although 
the detailed dependence on the shape of the object may 
be quite different. Finally, we note that the discussion 
presented here is dealing with a two-dimensional system. 
The analysis can be reproduced with similar results for 



APPENDIX 

Since we should recover equipartition in absence of dis- 
sipation, we expect that the kinetic energy of the motor 
is of the order of M(V 2 ) oc kT. It is therefore conve- 
nient to switch to the dimensionless quantity x of order 
1, defined as: 



(15) 



The Master equation is equivalent with the following set 
of coupled equations, describing the time evolution of the 
moments (x n ) — J x n P(x,t)dx: 




d t (x) 
d t (x 2 ) 



ddx 6 ) = 



= 2(xA 1 (x)) + (A 2 (x)) 
3(x 2 Ai(x)) +3(xA 2 (x)) 



(16) 
(17) 

(A 3 (x)) (18) 



with A n {x) the so-called jump moment, defined as: 



A n (x) 




u n W{^Jkt/Mx-u)&u. (19) 



Their explicit expression reads: 
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(20) 



The function $ [a, 6, z] is the Kummer function. This new set of equations are fully coupled and equally difficult to 
solve as the original Master equation. However, the equations decouple in the limit e = y/m/M — > 0. Indeed, using 
the following series expansions: 



, n 1 
1 + -,-,z 
2' 2' 

n + 3 3 



one finds: 



f + (2 + n)z 2 + -(8 + 6n + n 2 )z 4 + ... 
6 

1 + ~(3 + n)z 2 + i(15 + 8n + n 2 )z 4 + 
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(21) 
(22) 

(23) 



and: 



A 2 (x)^(l+r) 2 Sp 



2 'sin 2 e)e 2 + -x(sin 3 e)e 3 

7T 2 



(24) 



The expressions © and (flU)) for (V) and (T^ 2 ) are now moments (a;) and (a: 2 ) at the steady state up to order e 3 . 
readily obtained, by solving the equations for the first two 
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